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Abstract 

The ttNN form factor g-wNNiq 2 ) is obtained from a quenched lattice QCD 
calculation of the pseudoscalar form factor gp(q 2 ) of the proton with pion pole 
dominance. We find that g^NN^ 2 ) fitted with the monopole form agrees well 
with the Goldberger-Treiman relation and is much preferred over the dipole 
form. The monopole mass is determined to be 0.75 ± 0.14 GeV which shows 
that g-KNNiq 2 ) is rather soft. The extrapolated irN coupling constant g^NN = 
12.7±2.4 is quite consistent with the phenomenological values. We also compare 
g-wNN{q 2 ) with the axial form factor gA(q 2 ) to check the pion dominance in the 
induced pseudoscalar form factor hA(q 2 ) vis a vis chiral Ward identity. 



PACS numbers: 12.38.Gc, 14.20.Dh, 13.75.Gx, 13.60.-r 



The irNN form factor g^NN^ 2 ) is a fundamental quantity in low-energy pion- 
nucleon and nucleon-nucleon dynamics. Many dynamical issues like irN elastic and 
inelastic scattering, NN potential, three-body force (triton and 3 He binding ener- 
gies), pion photoproduction and electroproduction all depend on it. Similarly, the 
pseudoscalar form factor is important in testing low-energy theorems, the chiral Ward 
identity and the understanding of the explicit breaking of the chiral symmetry. Yet, 
compared with the electromagnetic form factors and the isovector axial form factor of 
the nucleon, the pseudoscalar form factor gp(q 2 ) and the nNN form factor g w NN(q 2 ) 
are poorly known either experimentally or theoretically. 

Notwithstanding decades of interest and numerous work, the shape and slope 
of gnNN{q 2 ) remain illusive and unsettled. Upon parametrizing g^NN^ 2 ) in the 
monopole form 

g-KNN\Q ) — gwNN-r2 2" I 1 ] 

A wNN 1 

with g n NN = g-wNNijni), the uncertainty in the parametrized monopole mass A n NN 
can be as large as a factor 2 or 3. For the sake of having a sufficiently strong tensor 
force to reproduce the asymptotic D- to S- wave ratio and the quadrupole moment 
in the deuteron, A nNN is shown to be greater than 1 GeV [[[]. Consequently, A wNN 
in the realistic NN potentials are typically fitted with large A^nn (e.g A^at ranges 
from 1.3 GeV to 2.3 - 2.5 GeV 0). On the other hand, arguments based on 
resolving the discrepancy of the Goldberger-Treiman relation [[| and the discrepancy 
between the ppn and pmr + couplings |J suggest a much softer g n NN(q 2 ) with A^nn 
around 0.8 GeV. Furthermore, hadronic models of baryons with meson clouds like 
the skyrmion typically have a rather soft form factor (i.e. A nNN ~ 0.6 GeV) due 
to its large pion cloud and such a small A nNN is needed for the high energy elastic 
pp scattering . 

In view of the large uncertainty in g^NN^ 2 ), it is high time to study it with 
a lattice QCD calculation. Since our recent calculations of the nucleon axial and 
electromagnetic form factors are within 10% of the experimental results |Q, 0, a 
prediction of g^NN^O 2 ) with a similar accuracy should be enough to adjudicate on 
the controversy over the 7rNN form factor. In this letter, we extend our lattice 
calculation to the proton pseudoscalar form factor for a range of light quark masses. 
g-KNNi.1 2 ) is obtained by considering the pion pole dominance in gp(q 2 ) when the 
latter is extrapolated to the quark mass which corresponds to the physical pion mass. 

In analogy with the study of the electromagnetic and axial form factors |5], || of 
the nucleon, we calculate the following two- and three-point functions for the proton 

G£(t,fi = E^^(0|T( X Q (x)x a (0)|0) (2) 

X 

Gf Pp (t fl p,t,q) = E e-^f + ^(0\T( X a (x f )P(x)^(0))\0), (3) 
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where x a is the proton interpolating field and P(x) is the mean-field improved isovec- 
tor pseudoscalar current for the Wilson fermion 

P{x) = ^L e m ^(x)i l5 ^(x). (4) 

Here, we have included the 2k/8k c (k c = 0.1568 is the critical k value for the chiral 
limit for our lattice at /3 — 6.0) and the e mqa ( m q a = ln(4« c /K — 3) is the quark 
mass) factors in the definition of the lattice current operator. These factors take into 
account the mean-field improvement and finite quark mass correction for the Wilson 
action jllj and have been shown to be an important improvement in the evaluation 



of the axial form factor in order to allow the perturbative lattice renormalization to 
work |Q. 

Phenomenologically, the pseudoscalar current matrix element is written as 

(ps\P(0)\p's') = gp(q 2 )u(p, s)i<y 5 u(p' , s'). (5) 

where gp(q 2 ) is the pseudoscalar form factor. It has been shown [§, [| that when 
tj — t and t » a, the lattice spacing, the combined ratios of three-point and two- 
point functions with different momentum transfers lead to desired form factors related 
to the probing currents. In the case of the pseudoscalar current in eq. (|J), the lattice 
pseudoscalar form factor gp(q 2 ) is determined from the following ratio 

T^G a / Pp (t f ,0,t,q)G^(t,6) q 3 

G%(t ft 0) Gf°(t, g) ~ * E q + m 9p{q } 1 ' 

where r = ( °^ jj V and m and E q are the proton mass and energy with momentum 
q respectively. 

Quark propagators have been generated on 24 quenched gauge configurations on 
a 16 3 x 24 lattice aX (3 = 6.0 to study the nucleon electromagnetic and axial form 
factors || |J|. We shall use the same propagators for the present calculation. Results 
are obtained for three light quarks with k = 0.154, 0.152, and 0.148. They correspond 
to quark masses m q of about 120, 200, and 370 MeV respectively. (The scale a -1 = 
1.74(10) GeV is set by fixing the nucleon mass to its physical value.) Results of gp(q 2 ) 
for the momentum transfers q 2 a 2 = n(2n/L) 2 (n = 1 to 4) are obtained from the 
plateaus of the ratio in eq. (|6|) as a function of t, the time slice of the current insertion, 
away from the sink and source of the nucleon interpolation fields |J. Since the ratio 
in eq. @ is proportional to q%, gp(q 2 ) at q 2 = can not be obtained directly. Rather, 
it will be obtained from extrapolation from the finite q 2 data as explained later. 

Plotted in fig. 1 are the lattice isovector pseudoscalar form factors gp(q 2 ) of the 
proton as a function of the quark mass in dimensionless unit m q a which takes into 
account the tadpole-improved definition for the quark mass. They include different 
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momentum transfers with q 2 from 1 to 4 times (2tt /La) 2 (N.B. q 2 = (E—mN) 2 —q 2 for 
the four-momentum transfer squared). The errors are obtained through the jackknife 
in this case. The extrapolation of gp to the quark mass m q a which corresponds to 
the physical pion mass is carried out with the correlated fit to a linear dependence on 
the quark mass m q a for k = 0.154, 0.152 and 0.148. The data covariance matrix is 
calculated with the single elimination jackknife error for gp which takes into account 
the correlation among the gauge configurations || [II]]- This fitting gives x 2 I^df = 
0.005,0.008,0.65, and 1.7 for q 2 a 2 form 1 to 4 (2vr/L) 2 . 

To extract the itNN form factor g-nNNio 2 )-, we take the pion pole dominance in 
the dispersion relation for gp(q 2 ) so that 

m z_ _ gz 

where G w = (0|P(0)|7r) can be obtained from the two-point function 

(0\Y,mt)P(0M0) t ^ a ^e~^ (8) 

Plotted in Fig. 2 is g^NNiQ 2 ) defined via eqs. (0) and There is a caveat to 

extracting g^NN (q 2 ) this way which we wish to point out. Strictly speaking eq. (0) is 
equivalent to PCAC where the physical pion field dominates and is thus valid for small 
q 2 . For q 2 as large as m 2 , with it' being the radially excited pion at 1.3 GeV, higher 
mass contribution to gp(q 2 ) may not be negligible. However, g^NN is expected to be 
an order of magnitude smaller than g„NN due to the the fact that a node in the internal 
qq wavefunction of it' will lead to cancellation in the vertex function. Therefore, we 
estimate that the pion pole dominance (eq. (0)) may have an error as large as 5 to 
10% at the highest q 2 we calculated. This is much smaller than the statistical error 
we have at the highest q 2 . This is also consistent with the estimate that PCAC and 
chiral perturbation is good to a scale of 4irf n . Keeping this in mind, we discuss the 
behavior of g-nNNi.1 2 )- We fitted it with both a monopole form (i.e. eq. flip) and a 
dipole form. We found that the monopole form with A tNN = 0.75 ± 0.14 GeV and a 
X 2 ' /Ndf — 0.13/2 is only slightly better than the dipole form with a dipole mass of 
1.32 ± 0.17 GeV and a x 2 /Ndf = 0.57/2, in so far as the x 2 is concerned. However, 
we can inject our knowledge at q 2 = 0. The Goldberger-Treiman (GT) relation which 
relates coupling constants at q 2 = 

m N g A (0) = APttjvjv(O) (9) 

predicts g vNN (0) = 12.66 ± 0.04 from the known g A (0) = 1.2573 ± 0.0028 and f n = 
93.15±0.11 MeV fll2| . Extrapolation of the monopole fit of g n NN(<l 2 ) gives g n NN(ty = 
12.2 ± 2.3 which agrees with the GT relation. Yet, the dipole fit giving gvjvjv(O) = 
10.8 ± 1.3 falls outside the prediction of the GT relation. Thus, we conclude that 
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the monopole form is much preferred over the dipole form. The monopole mass 
A-knn — 0.75 ± 0.14 GeV thus obtained is much smaller than those typically used in 
the NN potential, but agrees well with those based on the consideration of the GT 
relation [|J, the apparent discrepancy between g^ pp and g^+np ||, and the nucleon 
models like the skyrmion || . To salvage the nice fit of the NN scattering data and the 
deuteron properties based on a hard irNN form factor, attempts have been made to 
incorporate a soft g n NN(<l 2 ) either by appending a heavy pion at ~ 1.2 GeV [^3j or by 
including multi-meson exchanges JTJ] (e.g. Tip and ncr). Extrapolating g n NN(l 2 ) to 



q 2 = m 2 , we obtain g-^NN, the nN coupling constant, to be 12.7 ±2.4. This compares 
favorably with the empirical value of 13.40 ± 0.17 pfl and 13.13 ± 0.07 p[. The 4% 



change in g^NN^ 2 ) from q 2 = to m 2 indeed can account for the 4% discrepancy in 
the GT relation when the physical g n NN is used in eq. (|9|) instead of the (^ata^O) 0. 

Putting the chiral Ward identity d^A^ = 2m^i r y^T a /2^! with pion pole dominance 
or equivalently PCAC {d^A® = f^m 2 ^) between nucleon states, we find 

O I 1\ , 2r / 2\ 2 m lfn gnNN{q 2 ) / 1rA 

2m N g A {q ) + q h A {q ) = ~ 9 5 (10) 

m z _ gz 

In addition to PCAC, if one further assumes that the induced pseudoscalar form factor 
h,A{q 2 ) is dominated by the pion pole, i.e. /^a(<? 2 ) = 2f 7T g 7T N n (<? 2 )/( m ^ — <? 2 )> then 
g-KNN(<i 2 ) = (ttin/ fw)gA(q 2 )- In other words, g n NN(q 2 ) has the same q 2 dependence as 
5 , a(? 2 ) which has been frequently used in the literature [|T7], 0. As there is no a priori 
reason why g^NN^ 2 ) should have the same falloff as gA^Q 2 ) at all q 2 and, furthermore, 
chiral perturbation calculation (T^j at one loop suggests that they acquire different 
contributions, we compare g n NN(<l 2 ) from eq. (0) and gA^Q 2 ) obtained on the same 
set of gauge configurations |§ for the light quark cases. Both g^NN^ 2 ) and gA^Q 2 ), 
normalized at q 2 = 0, are plotted in Fig. 3 for k = 0.148,0.152,0.154, and 0.1567. 
The last n corresponds to the physical pion mass. We find that in all these light 
quark cases, there is a tendency for the normalized g n NN{<l 2 ) to lie lower/higher than 
the normalized gA{(f) at lower/higher — q 2 . This presumably reflects the preferred 
monopole vs dipole fit for the (yVTVA^g^and gA^q 2 )- Our data do not discern this well 
though. If this behavior is verified, it would imply that the induced pseudoscalar form 
factor hA^q 2 ) (not the pseudoscalar form factor gp(q 2 )) is not entirely dominated by 
the pion for higher — q 2 as it is at very low — q 2 (< 0.1 GeV 2 say). 

Lastly, from the chiral Ward identity (eq. (TTO])), we can obtain the induced pseu- 
doscalar form factor ^(g 2 ) from gA^Q 2 ) @ and g^NN^I 2 ) ■ The pion decay constant 
/„- needed in eq. (|T0D is calculated from the two-point functions (S^A 4 (t, x)P(0, 0)) 
and eq. (|8]). It is found to be 89.8 ± 4.5 MeV when the finite lattice renormaliza- 
tion is taken into account. We plot ^(g 2 ) in Fig. 4. Also plotted in the insert are 
experimental data obtained from pion electroproduction It turns out that the 
momentum transfer ranges of our lattice calculation and the available experiment do 
not overlap. We can not compare them directly. However, if we use the monopole 
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fit of g-xNNi.q 2 ) and the dipole fit of gA^q 2 ) [0, we find the extrapolation of hA^q 2 ) 
(solid line in Fig. 4) does agree with the experimental data at small — q 2 . We note 
the errors of the fit start to diverge as — q 2 — > due to the q 2 singularity in eq. ([10|). 
As a result we are not able to extrapolate to —q 2 = 0.88m 2 to compare with the the 
muon capture experiment. 

To conclude, we have calculated the isovector pseudoscalar form factor of the 
nucleon in a lattice QCD calculation for quark masses from about one to about two 
times that of the strange quark. From these we extracted g n NN(q 2 ) with the help of 
the pion pole dominance. The main results we gleaned are the following: 

1) Incorporating the Goldberger-Treiman relation at q 2 = 0, we find that g^NN^ 2 ) 
is much better described by a monopole than a dipole form. The monopole mass A^at 
= 0.75 ± 0.14 GeV is much smaller than commonly used in the NN potential. 

2) g-wNN = 12.7 ± 2.4 agrees well with the phenomenological values of 13.40 ± 



0.17 |T5| and 13.13 ± 0.07 fll6| . It is also consistent with the lattice calculation of 
14.8 ± 0.6 with staggered fermion [2~U|. 

3) The falloff of g n NN(q 2 ) is about the same as 5u(<? 2 ) at very small — q 2 (< 
0.3 GeV 2 ), but is likely to fall slower at higher —q 2 . This suggests that the induced 
pseudoscalar form factor tiA(q 2 ) is not entirely dominated by the pion pole at higher 
— q 2 . This point needs to be verified further with higher statistics study. 

4) From the chiral Ward identity and PCAC, we obtain h^iq 2 ) which can be 
checked experimentally in the future. 

The soft g n NN{q 2 ) form factor agrees with the predictions based on the discrepancy 
of the Goldberger-Treiman relation [|J and between the ppn and pmr + couplings []5J. 
This will have a large impact on the study of NN potential, the three-body force, and 
other processes which involve the ttN coupling. For future studies, it is essential to 
improve the calculation by expanding the volume in order to access smaller — q 2 and 
incorporating dynamical fermions effects. 

This work is partially supported by DOE Grant No. DE-FG05-84ER40154 and 
NSF Grants Nos. STI-9108764 and PHY-9203306. 
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Figure Captions 



Fig. 1 The lattice isovector pseudoscalar form factors at various — q 2 as obtained 
from eqs. (fj]) are plotted as a function of m q a, the quark mass in lattice unit, for the 
three light quark cases (Wilson k = 0.148,0.152, and 0.154). The top curve is for 
q 2 = (2ii/La) 2 , the rest are for q 2 from 2 to 4 times of (2n/La) 2 in descending order. 

Fig. 2 QnN n ((7^) at the quark mass which corresponds to the physical pion mass. 
The solid and the dashed curves represent the monopole and dipole fits with the 
respective monopole and dipole mass A. They give quite different extrapolations at 
q 2 = 0. 

Fig. 3 Comparison of gA^q 2 ) and g^NN^ 2 ) (both normalized to 1 at q 2 = 0) as a 
function of — q 2 for the light quark cases (k = 0.148, 0.152, 0.154, and 0.1567). 

Fig. 4 The induced pseudoscalar form factor /ia(<? 2 ) from eq. (|T0|) . The solid line 
is from the fits to <m(<? 2 ) and g n NN(q 2 )- Also plotted in the insert are data from the 
electroproduction of pion ||18|| . The typical size of the error bars for the solid line is 
indicated in the insert. 
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